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NOMENCLATURE 


Weibull distribution shape or scatter controlling 
parameter 

Weibull distribution scale parameter or characteristic 
life 

Weibull distribution scale parameter of time to first: 
failure 

Cumulative distribution of time to first failure 
Cumulative distribution of population 
test sample size 

number of failure observations in the test sample 

fleet size 

fleet reliability 

scatter factor for a specified reliability level 
time to failure variable 
time to first failure 

safe life for a specified reliability 
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1. INTRODUCTION 

Conventional procedures of reliability estimation are 
based on the assumption of arbitrary scatter factors of between 
2 and 4, which applied to an "estimated" life determined from 
development tests, are supposed to produce the safe and certi- 
fiable life. The probability of occurrence of individual lives 
shorter than the design life in the fleet of airplanes to 
which this estimate refers remains thus completely undefined 
but is presumably zero. No statements concerning the reli- 
ability level of the design can be made, since no confidence 
level can be attached to the selected scatter factor, nor can 
any statistical significance be attached to the "estimated" 
life. 

A basic improvement of the conventional procedures based 
on arbitrary scatter factors applied to an ill-defined "ex- 
pected life" was achieved through the application of the con- 
cept of the "expected time to first failure in a fleet" of 

structures subject to progressive fatigue damage, as shown 
1 * 

by Freudenthal . This research formed the basis of a further 

2 

investigation by Whittaker and Besuner on a reliability 
analysis plan for application to aluminum alloy structural 
fatigue performance which was compared to the current fixed- 
scatter-factor procedure for determining the safe life of a. 
structural detail. 


References at the end of report in numerical order. 
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The present study attempts to show the interrelation between 
the various parameters, such as fleet size, shape-parameter, 
level of reliability, in their influence on the scatter factor. 

The third asymtotic extremal distribution (Weibull) of 
smallest values was used for estimation of the certifiable life. 


2. MATHEMATICAL FORMULATION 

2.1 Estimation of Characteristic Life 


With the shape-parameter assumed to be known, one depends 
on the full- scale fatigue test to obtain the most likely value 
of the parer t population's central tendency or scale parameter. 

It is to be recognized that the occurrence of high-time outliers 
in this testing could lead to overly high safe-life estimator 
It can be argued that these high-time outliers are not part of 
the parent population and can therefore be censored. Fortunately, 
as has been observed in the data, large, complex, multicomponent 
specimens are not characterized by outliers to anywhere near the 

extent that simple monolithic specimens are. 

3 

As shown by Mann , the Maximum Likelihood Estimator (MLE. 
provides a best estimate for the two parameter Weibull distri- 
bution. 

. . . 2 3 

The point estimate of the Weibull characteristic life '3, 
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which considers the case of single-stage censoring at the 
failure time only. Equation (1) gives an unbiased uniformily 
minimum variance estimator for which the sampling distribution 
is known. For uncensored or complete sample, equation (1) 
reduces obviously to 

, n f 1/a 

3 = [-i ( Z Y®)] (2) 

n f i=l 1 

For this purpose it is also possible, if desired, to calculate 
' the lower bound interval estimate. The sampling distribution 

A 

of the Weibull 3 estimator is then given in terms of a prob- 
ability statement of confidence level, accounting for the 
finite sample size. 

2 . 2 Estimation of Certifiable Life 

If avoidance of a single catastrophic failure is a design 
requirement, the time to first failure will be the significant 
design criterion. 

From basic order statistical considerations- 1 ' 


1 " F Y 1 ( yi ) = U - F y (y)] 


N 


and 


a 


1 ~ F (y ) = exp [-N(q— ) ] 


1 - F (y ) = exp [ - 1 

Y 1 1 3 1 


(3) 


(4) 


3 




_1 

cx 

where = gN , which shows that for this distribution the 
function (1 - Fy^(y^)) is °f the same extremal type as (1 - F y 

with identical shape-parameter a but characteristic value 

-1 

reduced by (N a ) . Equation (4) can be solved for y , 






•)} 


1/a 


(5) 


Setting 1- F. (y, ) =R, the equation above can be rewritten 
1 1 




1/a 


(6) 


where Y R is the safe life for a specified fleet reliability. 
Equation (6) can also be called the "inverse reliability 
function." 

2.3 The Scatter Factor 

The scatter factor accounting fcr both the fleet exposure 
and specified fleet reliability R is simply 


s = £ 
b R Y 


(7) 


R 


as shown in figure 1. 

3. COMPUTATIONAL PROCEDURE 

In order to perform a parameter study a general computer 
program was developed and listed in Appendix A. The necessary 
input information is summarized in Table 1, 

The program can handle both censored and uncensored 
samples resulting from fatigue tests for the evaluation of 

A 

the MLE of 6. 


(y) ) r 
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Various shape parameters a can be specified according to 
the types of materials being used. Furthermore various 
fleet sizes can also be specified. As final input parameter 
the various levels of desired reliability can be specified. 

4. RESULTS OF DATA ANALYSIS 

As an example, item 56 of reference 2 was chosen at 
random which is an uncensored sample of size 3. The input 
information, is shown in table 1. The values for the shape 
parameters a vary from 1.0 up to 15.0, the fleet sizes from 
5 up to 5000. Levels of reliability range from 0.50 to 0.99. 
Tables 2 through 13 summarize the results in terms of scatter 
factors . 

5. DISCUSSION OF RESULTS 

The point estimate of the central tendency or characteristic 
life of the two-parameter Weibull distribution shows an incrsas- 
ing confidence in the higher life test results with increasing 
a (figure 2) . This tendency is quite obvious since a large 
value of a for the Weibull distribution corresponds to a small 
scatter in data since a(JtnY) = ir/a/6. 

The study also reveals the influence of fleet size, level 
of reliability as well as the Weibull shape parameter on the 
scatter factor. 

Figures 3 to 7 show the relationship of the reliability R 
to the scatter factor S R for the fleet sizes N = 5, 50, 250, 

1000 and 5000 respectively. Each plot shows a family of curves, 
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each curve corresponds to a certain a. One can note a con- 
siderable faster increase of scatter factors in the reliabil- 
ity interval of 0.95 and 0.99 than in the lower intervals. 

The results furthermore show a rapid increase in the 
values for scatter factors with increasing fleet size and level 
of reliability. This can be illustrcited by a simple example. 

For a metal with a predetermined a of 2 and a desired 

reliability level of 0.75 the scatter factor S n for a fleet 

of N = 5 will be 4.27; for a fleet of N = 250, S-, will already 

K 

be 29.50 and for N .= 5000, S D = 131.83. If a level of relia- 
bility of 0.99 is desired, the values for S R increase to 22.3, 
157.72 and 705.33 respectively. 

The scatter factors for metals with shape factors a < 2 
at higher levels of reliability (> 0.90) reach extremely high 
values especially for larger fleet sizes. See tables 2 to 5 
and figures 3 to 7 . 

Figures 8, 9 and 10 show a logarithmically linear rela- 
tionship between fleet sizes N and scatter factors S R for 
various levels of reliability (0.90, 0.75 and 0.50). For 
example: on the basis of a reliability of 0.90 a metal with 
a = 4 will yield higher values for scatter factors than a 
metal of a = 3 on the basis of a reliability of 0.5. But 
this will be the case only up to a fleet size of N = 200; if 
N > 200 the reverse is true. 
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Finally in figure 11 the relationship of the Weibull shape 
parameter and the scatter factor is demonstrated. The figure 
shows various families of curves. Their constant parameters 
are the fleet size N and the reliability R. It is interesting 
to note that for a given shape parameter a for example a fleet- 
size of N ~ 250 with a level of reliability of 0.50 results in 
lower scatter factors than a fleet size of N = 50 with a level 
of reliability of 0.90. The same relationship is true for both 
fleet sizes N = 250 and 1000, also for N = 1000 and 5000. 

Figure 11 also implies that the conventionally used scatter 

factors S R <_ 4 associated with an acceptable reliability level 

R = 0.9 would require a shape parameter a >_ 4 . 5 even for a 

fleet size N = 50 and a _> 6 . 0 for a fleet size N = 250 that is 

more representative of real conditions. Shape parameters of this 

magnitude imply a level of material-, production-, and mission- 

control that does not reflect realistic conditions of aircraft 

manufacture and utilization. Nor does a reliability level 

R = 0.90 reflect the implied expectation of the user. Thus, 

2 

under the assumption that a = 4 is a representative value for 
aluminum structures and that a reliability level R = 0.95 better 
reflects the expectation of the user, a scatter factor S R = 8.4 
(rather than S = 4.0) would define the "certifiable" life for 

K 

realistic fleet size of N = 250. This value of the scatter 
factor should therefore be used for long range transport planes 
of current design for which a = 4 is a representative value. 
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For advanced designs using titanium or high-strength steel 
alloys preliminary observations suggest that representative 
values of a are below a = 4', so that larger values of S_, 
for the specification of a certifiable life must be expected. 

6. CONCLUSIONS 

(a) The study shows clearly that a safe and economical 
design for metals with shape parameters a < 2 does not seem 
to be feasible if the design is based on the realistic con- 
cept of time to first failure, except for very small fleet 
sizes and low levels of reliability. 

(b) Scatter factors for establishment of a safe life 
used in procurement of fatigue sensitive aircraft are too 
narrow even for current designs. 
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Table 8 ; Scatter Factors for 
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Table 12: Scatter Factors for 
























Table 13: Scatter Factors 
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FIGURE 2 CHARACTERISTIC LIFE M 
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FIGURE 5 SCATTER FACTORS FOR FLEETSIZE N=250 



SCATTER FACTORS FOR FLEETSIZE N=1000 
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FIGURE 8 SCATTER FACTORS FOR RELIABILITY R=0.90 
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FIGURE 9 SCATTER FACTORS FOR RELIABILITY R=0.75 
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FIGURE 10 SCATTER FACTORS FOR.RE LIABILITY R=0.50 
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*** VARIABLES 


C SF = SCATTER FACTOR 

C. GAMA = C O N F 1 0 E N C E INTERVAL FOR FIR S T FAILURE 0 I S TR I BUT 1 O N 

C NG = NUMBER OF GAMMAS USED IN PROGRAM 

C FN = FLEET SI 2 E 

C NFN = NUMBER OF FLEETS I ZES 

C NS = SAMPLE S I ZE 

C NF = NUMBE 3 OF FAILURES ’ 

C ALFA = WEI3JLL SmAPE PARAMETER ! 

C Y = FATIGUE TEST SAMPLE 

C BETA = ML' . : OF WE I BULL SCALE PARAMETER ( POPULATION) UNBIAS ED 

C BETA1 = ' M I. E OF NEIbULL SCALE PARAMETER (FIRST FAILURE I UNBIASED 

C NA = NUMBER OF ALFAS ; 

C ZR1 = ESTIMATED LIFE FOR SPECIFIC CONFIDENCE INTERVAL (FIRS’ FAILURE ) 

C YF = FATIGUE FAILURE SAMPLE 

C YNF = NF FAILURE TIME 

C XI = TIME fO FIRST FAILURE __ 

C 

C *** DECLARAT 101 IS; ; ; 

C 

INTEGER FN ; 

DIMENSION ,U 1GCI ,B< ICC) 

D 1MENS 1 ON GAMA ( 20 ), ZR1 ( 200) . BETA 1( 2 CO I , BE T A t 2 00 ) 

DIMENSION FISH 50) , YF (20 ) , Y ( 20 ) . ALF A ( 60 I 
DIMENSION ;F(2Q0) . w(20l ; 

c 

C**-* INPUT - INFORMATION ; 

C 

REAP (5, 50 ' )) NS.NF,NA,NG.NFN 

500 FORMAT (5i:-,0) 

WRITE (6,5: Hi ; 

501 FORMAT ( 1 H ..,/////, 65X ,' I N P U T ',///) 

WRITE ( fc i 5 ( 1 2 ) NS , In E , f\ A , NG , NF N 

50 2 FORMAT (5X, ‘ SAMPLE SI At = ' , I 3, 3x , 'huMolR E A I L f 0= ' , 1 3 , ?X , ' AL F AS = ' , 

XI 3 . 3x , « GfaMnS= 1 , ! 3 , 3x, 1 in UMBER UE F L F I TS = 1 . I 3 , / / / ) 

READ (5, 50a) I Y< I ) , 1 = 1 ,NS) 

50 3 FORMAT (10. ; , ’ S A M P L E S 1 , / ) _ 

WRITE (6,503) 

WRITE (b,5i 'h ) ( Y( I ) . i = 1 , NS ) 

504 FORMAT (1CEb.1I 

_R£ ACi — Lli-5 If * ) — ( YE ( 1 ) , 1=1, NF ) . 

505 FORMAT 1//.1CX,' FAIL U R E S',/) 

WRITE (6,505) 

WRITE (6,5 04 ) ( YF ( I ) , I = 1,NE) 

READ (5,700) YNF 

700 FORMAT (F10.1) ^ 

READ 15,506) ( A 1 F A ( I ) , I = 1,NA) 

wRITE ( 6 , 520 I 

520 FORMAT (//.IQx.'A L E A S'. /I 

50b FORMA! (6F1C.2) ' 

WRITE (6,506) ( A L F A ( I ) , I = 1,NA) 

WRITE. (6,507) 

507 FORMAT (//.IQx.'F L E L- T S I Z E S',/) 

READ (5,501) (FIM(l), I = l,i\r,\) 
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31 


508 

FORMAT 

(1615) 


32 



wRlTE 

( 6,506) ( F N < I ) , 1 = 1, i.FN ) 


33 



WRITE ( 61 509) 


34 


5 09 

FORMAT 

( // / ,6X, 'GAMA ' , 7x , ' W , / / ) 


35 



DO 1 I 

= 1,NG 


"36 



READ (5,510 GAMA ( I ) ,-w ( I ) 

-\frv - 

37 


510 

FORMAT 

(2F10.3) 


38 



WRITE 

(6,511) I ,GAMA< I ) , *( 1 ) 

i 

39 


511 

FORMAT 

( 1 X , I 2 , F 7 . 3 , F 1 C . 3 ) 


4C 


1 

CONI INUE 



C 




TlV 


C 

T 5 # 

6 S T 1 MAT E ME I BULL SCALE PARAMETER 



C 






C 

C 

4(1) UN It NS OR ED SAMPLE* 


41 



00 2 j 

= 1 * N A 


42 



WRITE 

(6,516) ALF A (J) 


43 


518 

FORMAT 

MMl ,// / ,65 X, • ALFA= • , F6 . 2, / ) 


44 



IF (NS 

• N E % N F ) CO TO 4 


45 



SUMY = 

c.o 


46 



DO 3 L 

= 1 , N S 


47 



SUMY 

= SUMY + ( Y (L ) * * A L F A ( J ) ) 


48 


3 

CONTINUE ' 


49 



BETA(J) = (l./NF*SUMY I **( i. / ALFA < J ) ) 


5 0 



WE I T E 

(6,512) dETA(J) “ 


51 


512 

FORMAT 

( //,1CX, ' 6 E T A = ' ,F1G.1,//) 


52 



GO TO 

5 


53 

C 

4 

CONTINUE 



C 

c 


(2) Z ENSOP.EO SAMPLE 

1 


54 



SUMY = 

0 . 0 1 


55 



DO 6 I 

= 1 , NS 


56 



Sumy 

= SUMY + ( Y(LM*ALFA(J) ) 

' 

57 


6 

continue 


58 



6ETA(J) =1 i./NF*( SUMY t ( (NS-NF ) *( YNF**ALF A ( J ))).))**( 1 . /AL F A ( j ) ) 


59 



WRITE 

(6, ‘>191 BETA ( J ) 


60 


519 

FORMAT 

( // , 1 0 X , ' b E T A = 1 , F 1 0 . 1,//) 


fcl 


5 

CONTINUE 



~~<T 

c 

** * 

CALC.UL 

ATE Xl(GAMA) 


62 

c 


WP. 1 T E 

(6,515) (FN( I) , 1=1 , NFM 


6 3 


51~5 

FORMAT 

( / / ,3x, » GAMA' ,7 a, 'Ms* , I 2 , 5x , 1 ri= ',13,5a, 'K = ' , I 3 , 4 a , ' i\ = '., I 4 , 





X4X , ' N= 

' , 1 4 , 4 X , 1 N= ' , 1 4 , 4 X , •;>!=' , I , 4 < ,’• f\' = ' , I 4 , H X , ' N= ' , 1 -t , -tX , ' :i= ' , 





X I 4 | / / / ) 


64 



DO 8 N 

= 1 , N G 


65 



DO 7 I 

= ) ,NFN 


66 



BE T A 1 ( I ) - iiETA(J) * ' 1./ ( F.N( I 1 * * ( 1 . /ALF4( J ) ) ) 


67 



AIK) = 

(1-0 - GAMA ( K ) ) 


6 8 



B ( J ) = 

1.0 / ALFA(J) 


69 



ZR 1 ( I ) 

= B ET A 1 ( I M ( ( Ac S ( A LOG ( A ( K ))))*'( B ( J )) ) 


70 



SF ( I I 

= 3ETA( JI/ZRK I ) 


71 


7 

CONTINUE 


72 



WRITE 

(6,5131 GAMA (K) , ( SF (I) ,.I = 1 , N FM 


73 


TTI 

F OHM AT 

( F 3 . 3 , 1 C ( F 1 6 . 3 ) ) 


74 


8 

CONTINUE " ' 


75 


2 

CONTINUE 


76 



RETURN 
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END 
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